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We derive a master equation for a mirror interacting 
with the vacuum field via radiation pressure. The dynami- 
cal Casimir effect leads to decoherence of a 'Schrodinger cat' 
state in a time scale that depends on the degree of 'macro- 
scopicity' of the state components, and which may be much 
shorter than the relaxation time scale. Coherent states are 
selected by the interaction as pointer states. 

Within the framework of Quantum Mechanics, a 
closed system may be found in any quantum state of the 
Hilbert space. As pointed out by Schrodinger [|l[, this 
is in apparent contradiction with the classical behavior 
of macroscopic systems. However, macroscopic systems 
are seldom isolated, and the interaction with the environ- 
ment engenders the decay of most states into a statistical 
mixture of 'pointer states,' which are linked to classical 
properties of the system H. Coherent superpositions of 
pointer states decohere into a statistical mixture in a time 
scale which is usually of the order of the damping time 
divided by some parameter representing the degree of 
'classicality' of the states. The decoherence time scale 
for a microwave field in a high-Q superconducting cavity 
was recently measured 1^] to be in agreement with such 
prediction Q. 

Several different heuristic models for the coupling 
with the environment have been considered |5[|. In this 
letter, we show that the coupling with the quantum vac- 
uum field via radiation pressure provides a more funda- 
mental, ah initio model for decoherence. The Casimir 
effect for moving boundaries has attracted a lot of in- 
terest recently [p|. The vacuum radiation pressure force 
dissipates the mechanical energy of an oscillating mirror, 
and the associated photon emission effect could in princi- 
ple be measured experimentally M . Usually, one assumes 
that the mirror follows a prescribed trajectory, thus ne- 
glecting the recoil effect. However, here we want to focus 
on the mirror as a dynamical quantum system, hence 
the need to take the full mirror-plus-field dynamics into 
account. Jaekel and Reynaud treated this problem by 
using linear response theory B , in order to calculate the 
fluctuations of the position of a dispersive mirror driven 
by the vacuum radiation pressure. Mass corrections were 
also obtained in Refs. y] and flOfl . 

In this letter, we consider a nonrelativistic partially- 
reflecting mirror of mass M (position q and momentum p) 
in the harmonic potential of frequency wq, and under the 
action of vacuum radiation pressure. We take a scalar 
field in 1 -|- 1 dimensions, and neglect third and higher 
order terms in v/c, where v is the mirror's velocity (we 



set c = 1). We start from the Hamiltonian formalism 
developed in Ref . flo] . The Hamiltonian is given hy H = 
Hm + Hf + ^int, where 
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is the Hamiltonian for the scalar field (H = 9(0 is its 
momentum canonically conjugate) under the boundary 
condition corresponding to a partially-reflecting mirror 
at rest at a; = 0, where the coupling constant fl also 
plays the role of a transparency frequency, the frequency- 
dependent reflection amplitude being R^lo) = —in/{io + 
iQ) Ijl^l- Since the emitted photons have frequencies 
smaller than ujq, the perfectly- reflecting limit corresponds 
to ujQ <^ Q. We allow in principle for arbitrary values of 
LUo/n, but assume from the start that Tiujo/M, which is of 
the order of the recoil velocity of the mirror, is very small. 
The interaction Hamiltonian iJint describes, on one hand, 
the modification of the boundary condition for the field 
due to the motion of the mirror, and, on the other hand, 
the modification of the mirror's motion engendered by 
the radiation pressure force. The first effect leads to the 
emission of photon pairs out of the vacuum state (dynam- 
ical Casimir effect), whereas the second leads to dissipa- 
tion and decoherence of the mirror's motion, as shown 
below. To second order in v/c, we have 



H\, 



^^ + ^-V(0)4, (3) 



where V = — J dxdx(t>dt4> is the field momentum oper- 
ator. In the r.-h.-s. of (^, the first term is the most 
important, yielding the effects of dissipation and deco- 
herence. The second term does not depend on the mir- 
ror's variables, and hence will be of no relevance here, 
whereas the third term, being already of second order in 
v/c, is taken only to first order in perturbation theory. 
As discussed in Ref. 



|10 



it provides a contribution to 
the mirror's mass shift. 

We derive a master equation for the reduced density 
matrix of the mirror p{t) by assuming that at i = 
the mirror and the field are not correlated, so that the 
density matrix of the combined system p is written as 
(5(0) — p{0) ® pf, where pF is the density matrix of the 
field alone. Then we compute p{t) up to second order in 



the perturbation Hamiltonian iJmt- Note that the smaU 
perturbation parameter is the mirror's velocity ti/c, and 
not the couphng constant fl, which is incorporated in the 
field Hamiltonian Hp. In fact, one may diagonalize Hm + 
Hp for arbitrary values of f2, writing field eigenfunctions 
containing reflection and transmission coefficients so as to 
satisfy the boundary conditions at a; = |l0[ . Finally, the 
master equation for p{t) is obtained by tracing p{t) over 
the field variables, taking the field to be in the vacuum 
state. We find 



ihp = [Hm - 



AM{t) p2 
M 2M 



,p]~r{t)[p,{q,p}] 



-^D,{t)[p,[p,p]]-^D2{t)[p,[q,p]]. 



(4) 



The mass shift in (|) is given by AM{t) = AMi + 
AM2(i), where the cut-off dependent AMi =< fi<?!)2(0) > 
is the only (first order) contribution of the p^ term in 
Eq. (0). It was derived earlier by different methods in B] 
and |10| . Except for AAfi, the terms in (^) come from 
second order perturbation theory. The corresponding co- 
efficients are calculated from vacuum correlation func- 
tions of the momentum operator. The mass shift AA/2(t) 
and the damping coefficient r(i) are obtained from the 
anti-symmetric correlation function ^(i) = {['P{t),'P{0)]), 
which is connected to the susceptibility describing how 
the field momentum is affected by the motion of the mir- 
ror (and the corresponding modification of the boundary 
conditions). In fact, we show below that T{t) is closely 
connected to the photon emission effect and the associ- 
ated radiation reaction force that damps the motion so 
as to enforce energy conservation. The diffusion coeffi- 
cients Di (t) and D2 (t) are obtained from the symmetric 
correlation function cr(t) = {{'P(t),V{0)}) -2{'P)'^, which 
represents the vacuum fluctuations. 

Since V is quadratic in the held operators, the correla- 
tion functions are obtained from the two-photon matrix 
elements (017^(^)1^1, W2), which are calculated by using 
the normal mode expansion for the held operator. The 
spectral density S(w) is defined as the Fourier trans- 
form of S.(t). For a; > 0, '^(j-o) results from the contri- 
bution of two-photon states with uji + uj2 — uj. We find 
E{uj) = {2/TT)1\^nC{uj/n) with C(u) = ln(l + u^)/{2u) + 
(arctanu)/u^ — 1/w, whereas the Fourier transform of 
a{t) is e(a;)S(a;) (e(w) is the sign function). The trans- 
parency frequency fJ sets a frequency scale for the be- 
havior of S(a;). Thus, for a; <C f^ the spectral density is 
linear ('ohmic' environment), whereas for high frequen- 
cies it goes to zero as C(u) ~ \n.{u)/u, due to the mirror's 
transparency at frequencies w ^ 17. We find 
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The function sin [{ujq ^ uj)t] /{loq — oj) in equations (H) 
and (M) has a peak of width 27r/i at oj = wq- For large 
times, ilt 3> 1, the spectral density is approximately con- 
stant over the width of this peak, and then may be taken 
out of the integral, yielding 



r = ?^C(Wf^) - '"° 



27rM 



127rM' 



(9) 



the last approximation being valid in the perfectly- 
reflecting limit. If we also assume that uji^t ^ 1, Eq. (M) 
yields Di = fiT /{Muji^). Accordingly, for large times the 
damping and diffusion coefficients have constant values 
that result from the contribution of two-photon states 
\ijji,ijj2) such that LO1+UJ2 = uj = wq. This is precisely the 
condition satisfied by the photon pairs generated in the 
dynamical Casimir effect [|6|. In fact, the damping rate 
r as given by Eq. (||) is directly connected to the dissi- 
pative force on the moving mirror F = fix /(67r) |ll| ] 
(for simplicity we consider the perfectly reflecting limit). 
Indeed, the equation of motion for the average position 
then reads x = —uJ^x -\- hx /(GttA/), whose solution in 
the limit Tlloo/M <C 1 decays as exp[— ?iijjQi/(127rM)] in 
agreement with Eq. (||). 

The asymptotic values of the dispersive terms AM2 [t) 
and D2{t) do not originate, on the other hand, from the 
neighborhood of w = wq- In the perfectly-reflecting limit, 
we neglect wg in the denominator in Eq. (pi), and, when 
fit ^ 1, replace the sine squared by one- half. Integra- 
tion of the resulting expression over the whole frequency 
interval yields AAf2 ~ T%VI/{2'k). Accordingly, for large 
times we find the same mass correction obtained in llQ] 
from stationary perturbation theory. 

From these results, we may address two fundamental 
issues: (i) find out the pointer states; (ii) estimate the 
decoherence time scale. In the context considered here, 
pointer states are the most robust elements of the Hilbert 
space with respect to the motional interaction with the 
vacuum field. A simple test was proposed in Ref. [ p2[ , 
based on the idea that for pointer states the rate of in- 
formation loss is minimum. Such rate is measured with 
the help of the linear entropy s = 1 — Trp^ (s = for 
a pure state and greater than zero for a mixture). We 
calculate the rate of entropy increase from the master 
equation (^1, assuming that the initial state is pure: 



i(0 = 2r(t)(.(t) - 1) + l^(Ap)^ + ^, 



'li,p 



Wo ■ 



(10) 



where (Ap)^ is the momentum dispersion and <jq^p = 
({g,p}) — 2{p){q) (with all operators evaluated at the 
same time t). The first term in Eq. (|l^ leads to a de- 
crease of entropy (hence damping tries to localize the 
state competing against diffusion) which does not depend 
on the initial state. Thus, it is not relevant for the de- 
termination of the pointer states, and will be left out of 
our discussion. 
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FIG. 1. Difltusion and damping coefficients as a function of 
time in the perfect-reflector limit uio/fl = 10"* <g 1. Here 
Di = Tt^uio/Vl-KM^ and F = hujQ/12-KM are the asymptotic 
limits of D\(t) and F(f). The insets show the behavior of 
these two time-dependent coefficients for short times. 

We first consider the effect of the last two terms in 
Eq. ([lOl ) in the perfectly-reflecting limit. In Fig. 1 we 
plot the diffusion and damping coefflcients as functions 
of ijjQt for cjo/f^ = 10^^. Di{t) develops an initial jolt 
for times of the order of 51^^ and then decreases to the 
asymptotic value (-Di)pcrf = ^ ujo/{12t:1\''P) for t ^ l/^o- 
If we integrate Eq. (|o|) over many periods of oscillation, 
from t = to t = T = riZn Jloq, the contribution to the 
entropy of the initial jolt is negligible, allowing us to re- 
place the diffusion coefficients by their constant asymp- 
totic values. When computing (Tq,p{t) and (Ap)^(i) in 
Eq. (px|), we take the free evolution (corresponding to 
the harmonic oscillator Hamiltonian Hm) of the mirror's 
operators q and p (weak coupling approximation). We 
get 



s{T) = 2T^ [{Ap)l 



{MuJo)^Aq)l] 



(11) 



where (Ap)q and (Ag)Q represent the dispersions for the 
initial state. Note that D2{t) does not contribute to 
the time-averaged entropy production. The minimum 
s{T) given the constraint AqAp > h/2 is for Aq^ = 
h/{2MuJo), Ap2 = Mhuja/2. Thus, as in the problem 
of QBM with interaction Hamiltonian linear in the posi- 
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p2| , the pointer basis consists of coherent 
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FIG. 2. Diffusion and damping coefficients as a function of 
time in the high-transmission limit ujq/Q, = 10 ^1. Here 
Di = fi^^'^\n{ujo/Q.)/2-KM'^LJo and F = n^'^\n{ujo/Q.)/2-KM 
axe the asymptotic limits of D\{t) and F(f). 

The opposite limit o^o 3> ^i corresponds to domi- 
nant frequencies of the environment slow with respect 
to the mirror's own time scale. However, since the spec- 
tral density S(w) decays too slowly for u ^ Vt, field fre- 
quencies of the order of wq provide a significant contri- 
bution even in this limit. As a consequence, the vac- 
uum field does not behave as an adiabatic environment 
in the sense of Ref. Q. In Fig. 2, we plot the dif- 
fusion and damping coefficients as functions of uj^t for 
(jJq/VL = 10''. They oscillate around their asymptotic val- 
ues with (angular) frequency loq and with an amplitude 
of oscillation that decays in a time t ^ 1/fl ||l4|. The 
oscillatory terms do not contribute to the entropy in- 
crease when we average over many oscillations. Hence 
Eq. (^l|) also holds in this case, although the rate of 
entropy increase is much smaller than in the perfect- 
reflecting limit, since the asymptotic limit of Di{t) is 

now Di = 6(f7/wo)^ln(wo/^)(£'l)perf < (-Dl)pcrf- 

In order to estimate the decoherence time scale, we 
take, at i = 0, the 'Schrddinger cat' state \tp) = (|a) -I- 
I — a))/\/2, with a = iPo/^/2MhuJo- The corresponding 
Wigner function is 



W = W„ 



1 



exp 



2pHAq) 



21 



2(A9)2 



cos(^), 
(12) 



with Aq — y^'h/{2Mujo), and where Wm corresponds to 
the statistical mixture /?,„ — (l/2)(|Q:)(a| -I- | — a){—a\). 
In phase space, Wm has two peaks along the momen- 
tum axis (at ±Pq). The second term in Eq. (n2h origi- 
nates from the interference between the two state com- 



ponents, and hence represents the coherence of the cat 
state. Since it oscillates along the q axis in phase space, 
diffusion in position will damp the coherence at a max- 
imum rate given by —Did^W/W, with, from Eq. (12), 



d'^W 



Di{2Po/h)'^W. After averaging the decoher- 
ence rate over a period of oscillation p5[, we find that 
the decoherence time scale td is 
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To clarify the connection between decoherence and 
the dynamical Casimir effect, we present a second deriva- 
tion of Eq. (O), based on the concept of entanglement 
between mirror and field on account of the generation of 
photon pairs. At t = 0, the quantum state \'^) of the 
complete mirror-plus- field system is |^)o = IV") ® |0) In- 
stead of tracing over the field operators, we follow the 
evolution of the field state (in the interaction picture) 
to find |^)t = (|a) (g> |(/)+)t + \~a)(g) |(/)~)t)/\/2, where 
\(f> )t is computed from first-order perturbation theory 
assuming a classical prescribed motion: 



\^' 



where 



Bit)\0)± 



1 



duJi I duj2b{uJi,U2;t)\uJi,UJ2), 

(14) 



b{uji,uj2;t) = -{uji,uj2\'P\0) 



dt'e^("'i+"^)*'g(i') 



is the two-photon amplitude corresponding to the mir- 
ror's velocity q{t) = —iy^2Uujo/M a cos{uJot) associated 
to the state |a), whereas |i?(t)p is determined by the nor- 
malization condition (0 |0 ) = 1. Since the amplitude 
is proportional to the velocity, it has an opposite sign 
when associated to \ — a), as shown in Eq. (|lj) . When 
ujot 3> 1, the two-photon probabilities are proportional 
to the time t, and related to the relaxation rate F. Then, 
from Eq. (jl|) we derive p{t)-p,n = il-t/td){p{0)~ Pm), 
with td given by ([l3|). 

According to Eq. ([13|), decoherence is faster than en- 
ergy dissipation by a factor that represents the degree 
of 'macroscopicity' of the coherent states. In fact, |ap 
is twice the ratio between the energy of the coherent 
state and the zero-point energy of the harmonic oscil- 
lator. Therefore, Eq. ( [l3| ) provides an additional illus- 
tration of the meaning of the limit |a| 3> 1 as the classi- 
cal limit of the quantum harmonic oscillator. Moreover, 
Eq. (Ol) also shows that the decoherence rate increases 
with the distance beween the two coherent components 
in phase space. We have confirmed the role of coher- 
ent states in the understanding of the classical limit by 
showing that they are the pointer states. Remarkably, 
classical behavior is obtained from the mere inclusion of 
an unavoidable, intrinsically quantum effect, the radia- 
tion pressure coupling with the quantum vacuum field. 
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